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We study theoretically a binary system in which an attrac- 
tion of unlike particles is combined with a type-independent 
soft core repulsion. The possible experimental implementa- 
tion of the system is a mixture of DNA-covered colloids, in 
which both the repulsion and the attraction may be induced 
by DNA solution. The system is shown to exhibit surprisingly 
diverse and unusual morphologies. Among them are the dia- 
mond lattice and the membrane phase with in-plane square 
order, a striking example of spontaneous compactiGcation. 
PACS numbers: 82.70.Dd, 81.16.Dn, 87.14. Gg, 
42.70.Qs 

Self-assembly in colloidal systems has attracted a lot 
of interest as an experimental tool for study of crystal- 
lization and glassiness Q. It also has a great potential 
for the submicron technology, especially for fabrication of 
photonic band gap materials Q . The progress in these di- 
rections is considerably limited by the lack of diversity of 
the crystalline morphologies achievable by self-assembly. 
Typically, a monodisperse colloidal system crystallizes 
into a close-packed structure (FCC, or another stacking 
of hexagonal layers) , (|] . 

In this paper, we propose a system which combines 
relatively simple interactions with a rich and unex- 
pected phase behavior. It is inspired by recent experi- 
mental demonstration of DNA-assisted self-assembly of 
nanoparticles . The key elements of that scheme are 
submicron spheres (e.g. golden, silica, or other) cov- 
ered with short single-stranded DNA "markers". The 
marker sequence determines the particle type (there may 
be many markers per particle, but their sequences must 
be the same). One can now introduce type-dependent in- 
teractions between the particles by adding "linker" DNA 
molecules, whose ends are complementary to the corre- 
sponding markers. These interactions are very selective, 
reversible and tunable. 

We start our discussion with a generic model in which 
the physical origin and details of the inter-particle in- 
teractions are largely ignored. Let us consider a binary 
system of spheres, of the same diameter d, repelling each 
other with a soft-core potential, U{r). In addition, the 
unlike particles may stick to each other with binding en- 
ergy —E. We will study the phase behavior of the system 
for various values of two parameters: the aspect ratio, 
d/^ (^ is the range of potential U{r)), and the relative 
strength of the attraction, E/Uq (here Uq = U{d)). Later 
in the paper, we discuss how this system can be imple- 
mented experimentally. Our specific proposal is to use 
DNA to introduce type-dependent attraction, and poly- 
meric (possibly, also DNA) brush to induce the repulsive 



potential. 

The non-trivial phase behavior of the discussed sys- 
tem is a result of interplay between the adhesive energy 
and the soft-core repulsion. We will consider only struc- 
tures with 1 : 1 composition. Let be a distance to the 
fc-th nearest neighbor of a particle in a given structure 
(ri = d), let Zfe be the average number of such neighbors, 
and Z be the average number of cohesive contacts per 
particle (coordination number). If the entropic effects 
are neglected, the energy per particle (i.e. the average 
chemical potential of A and B particles) is 

li=\ (^-ZE + ZfcC/(rfc) j ~ i i-ZE + Z2U{r2)) 

(1) 

Here, we have neglected the contribution from the par- 
ticles beyond the second nearest neighbors, which is a 
reasonable approximation for d/^ 2> 1. In order to as- 
sure its validity, we perform a posteriori check of the 
effect of the higher order corrections on our results. It is 
straightforward to use the above equation to identify the 
phase boundary between two different structures. Since 
the chemical potential should be continuous at the transi- 
tion , one can express the critical value of adhesive energy 
E in terms of the geometrical parameters (Z , Z2, and r2) 
of the two phases 

Z^U{r'2) - Z2U{r2) 

z^ 

The task of identifying all plausible morphologies of 
our system is clearly more challenging than comparing 
them energetically. There is hardly any systematic way 
of doing this, which would go beyond an educated guess. 
Nevertheless, one can considerably limit the search by 
making a number of assumptions based on general prin- 
ciples j6|: (i) In order to avoid direct contacts of the 
same-type particles, the structures should be bipartal, 
i.e. they should consist of two sub-lattices correspond- 
ing to the two types of particles, so that all the nearest 
neighbors were of opposite types, (ii) All the nearest 
neighbors should have the same bond length. This is 
needed to take advantage of the cohesive energy, as long 
as we model the particles as rigid sticky spheres, (iii) 
The structure is likely to possess a high symmetry: we 
consider only crystalline morphologies, with all equiva- 
lent sites. The violations of this criterium are fairly rare 
(e.g. quasicrystals). 
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In the table below, the candidate phases satisfying the 
above principles are classified accordingly to their coor- 
dination numbers and the space dimensionality. 



z 


2D 


3D 


3 


Honeycomb 


Gyroid 


4 


Square (SQ) 


Diamond (D) 


5 




Honeycomb Stacking (HS) 


6 




Simple Cubic (SC) 


8 




Body-Centered Cubic (BCC) 



Note that in the limit of large aspect ratio, among var- 
ious phases with the same Z, the one with the largest 
is energetically preferred, independently on the number 
of the second nearest neighbors (Z2). Likewise, since 
goes down with the increase of coordination number, Z, 
the transition value of given by Eq. (||), asymptoti- 
cally reaches value Z2U{r2), determined by the parame- 
ters {Z2, 7*2) of the higher-Z phase. In other words, the 
higher Z, the higher is the parameter E{Z) of the tran- 
sition between Z- and Z + 1-coordinated phases. One 
can conclude that the sequence of the phases in the large 
aspect ratio regime is generic and should be indepen- 
dent of the particular choice of the repulsive potential: 
BCC(Z = 8)-SC(Z = 6)-Honeycomb Stacking {Z = 5)- 
Diamond(Z = 4). The possibility of the self-assembly of 
the diamond lattice is especially exciting, because of its 
potential as a photonic bad gap structure 1^. 

The energetic difference between the two Z = 3 phases 
in the above table is only due to interactions of the third- 
nearest neighbors, which makes it irrelevant for a realis- 
tic situation. Unless the system is confined to 2D, this 
suggests that the observed morphology with Z = 3 will 
typically be a disordered one. There is no true phase 
transitions for (Z < 3) because of the low dimensionality 
of the dominant structures. 

As the aspect ratio D/^ \s being decreased, one can 
expect two types of events: (i) "squeezing out" of cer- 
tain phases by their neighbors; (ii) replacingone struc- 
ture with another without changing Z. As one can see 
from Figure 0(a), the phase diagram obtained for expo- 
nential potential, U{r) = C/oexp(— (r — d)/S) provides 
examples of the both kinds. One of them is particularly 
striking: at certain aspect ratio, the system undergoes 
a transition from 3D diamond lattice to 2D membrane 
with in-plane square order. Even though this order is 
known to be strongly affected by long-range fluctuations 

, the corresponding Landau-Pierls effect does not have 
any divergent contributions to the chemical potential of 
the 2D phase. One may refer to D — SQ transition as 
spontaneous compactification. Although the found 2D 
phase is somewhat similar to lipid membranes, it is built 
by particles with isotropic effective interactions (unlike 
lipids). The same holds for the diamond lattice: unlike 
the diamond structures in nature, in our case it is due to 
the competition of relatively simple isotropic potentials! 
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FIG. 1. Calculated phase diagram of the system for expo- 
nential (a) and Gaussian (b) forms of the repulsive potential 
U{r). 

As it was argued, the major features of the obtained 
phase diagram are fairly independent of the choice of the 
form of the inter-particle potential. To check this, we 
have calculated the phase diagram for two types of in- 
teractions, exponential and Gaussian, U{r) — exp(— (r — 
(i)^/2^^) (see Figure The major difference is that the 
square lattice (SQ) completely disappears in the Gaus- 
sian case. In addition the diamond lattice (D) signif- 
icantly expands at the expense of SQ and HS phases. 
This trend appears to be quite general: the balance be- 
tween the two competing Z = 4 phases, D and SQ, shifts 
towards Diamond for potentials with a super-exponential 
decay, while the region of stability of SQ expands for 
sub-exponential U (r) (such as power laws, stretched ex- 
ponentials, or Yukawa interaction). 

We now proceed with the discussion of a plausible 
experimental implementation of the proposed system. 
As we have already mentioned, type-dependent "DNA 
bridging" of colloidal particles is an appealing way to 
introduce the AB attraction. As to the repulsive inter- 
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actions, there are several candidates in the colloidal sci- 
ence. Unfortunately, using of electrostatic repulsion in 
our case is problematic: in order for it to be relevant, the 
salt concentration needs to be much lower than at the 
physiological conditions. As a result, there would be a 
high electrostatic barrier for DNA duplex formation. 

An alternative way to introduce the soft-core repulsion 
is through steric interactions of polymer-covered parti- 
cles. Here we focus on a particular scenario in which the 
repulsion is due to DNA molecules with only one "sticky 
end" (complementary to either A or B markers). Un- 
like linker DNA, these one-arm molecules do not result 
in bridging and play a role of a buffer. We will assume 
that both the buffer-DNA and the linkers are double 
stranded, with the exception of the short terminal seg- 
ments. The strength of the interaction of these sticky 
ends with the complementary markers can be character- 
ized by DNA concentration cq, at which the condensa- 
tion would occurs (i.e. when the chemical potential of an 
adsorbed chain would become equal to that in the solu- 
tion). If the actual concentration of the buffer-DNA, c, 
is much lower than cq, the number of the adsorbed chains 
per particle is iV = A/'maxc/c^^^. Here N^nay^ is the total 
number of markers per particle. For the sake of simplic- 
ity, we assume that iV™'*'^ and N are the same for A and 
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B particles, i.e. ca/cb 

Note that we have totally neglected the excluded vol- 
ume interactions between the DNA molecules, which 
might seriously decrease the coverage compared to our 
result. The reason for this is twofold. First, the optimal 
coverage (discussed below) corresponds to a moderate 
overlap between the adsorbed chains. Second, the high 
rigidity of a double-strand DNA molecule results in the 
excluded volume effects being negligible within a single 
chain as long as it remains shorter than several thou- 
sands persistence lengths, Ip ~ 50 nm. Thus, the ad- 
sorbed DNA molecules can be treated as phantom Gaus- 
sian chains, which only interact with hard surfaces of the 
particles, but not with each other. Since the adsorption is 
reversible, the confined chains are being "squeezed out" 
when the gap between two particles, r — do becomes com- 
parable to the gyration radius of the DNA chain, Rg. The 
corresponding energetic penalty can be calculated in the 
spirit of the Deryagin approximation |0] (from now on, 
we distinguish bare particle diameter, do, and effective 
one, d): 



U{r) ~ 2NkT^ 

"0 J(r-do)/Rg 



1 — exp 



kT I 



dA 

(3) 



Here M^(A) is the free energy penalty for the confinement 
of a Gaussian polymeric chain between two walls at sepa- 
ration i?gA. It can be obtained by using the Schrodinger- 
like description of the ideal polymer p|,p|l: 



Ty(A) UVi^ 

\ Ji=0 



1 (Ti{l + 2ny 



6 V A 



(4) 



The resulting repulsive potential is shown on Figure ^ 




(r-d„)/Rg 

FIG. 2. Repulsive, U(r), and attractive , UAB(r) potentials 
induced by DNA-particle interactions. Solid lines correspond 
to R'g/Rg — 1. Note the barrier in the attractive potential for 
R'g/Rg = 0.5 (dashed line). 

The attractive potential induced by the linker-DNA 
can be calculated in a very similar manner. An impor- 
tant differences are that one has to take into account the 
elastic energy of a stretched linker, Wei{A) = 3kTA^/2, 
and that the condensation concentration of the chains 
with two sticky ends is given by — c'^c'g R'g. Note 
that in a general case the gyration radius of a linker- 
DNA, R'g, may be different from Rg. 



UAB{r)^U{r) 



2Nl^^kTcAB (R 



^AB 



(5) 



{r-d„)/R, 



exp 



WiA)+W^i(A) 



kT 



dA 



The relative strength of the attraction and the repul- 
sion is controlled by the ratio of the concentrations of the 
buffer-DNA and linkers: 



E 



CAB 

Cab 



- 1 



(6) 



Here Cab is the linker concentration at which E van- 



ishes. For the case of R'g ^ Rg, 



Cab 



o„(0)JO)o5 
jc^ rig 



(7) 



The absolute value of the energy scale in the problem 
is also controllable by concentration: according to Eq. 
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U„ ~ NkT^ 
a 



(8) 



This scale should be of order of several kT to ensure 
that the phase diagram is not affected by the thermal 
fluctuations, and that the escape time from the potential 
well of depth E = — min ([/^^(r)) is not too long. By 
comparing Eqs. (0) and (g), we conclude that Cab — 
CACBRy(P. 

As one might expect, the tail of repulsive potential 
U{r) is well described by a Gaussian with the character- 
istic length scale ^ ~ i?„: 



U{t) 



exp 



Rn 



(9) 



Here d is the effective diameter which is determined by 
the position of the minimum of UAB{f)i and the bias 
/3 ~ 0.15 + [d — do)/Rg. For the case R'g/Rg = 1, shown 
on Figure H, /3 ~ 0.6. Thus, the corresponding phase dia- 
gram should be somewhere halfway between the Gaussian 
and exponential ones. By increasing the ratio R'g/Rg one 
can move the system more towards exponential regime, 
because the position of the minimum changes roughly 
linearly with the radius of the linker DNA. However, 
the dynamic range of R'g/Rg is rather limited: the long 
linkers would result in the additional attraction beyond 
the nearest neighbors, which would violate our initial as- 
sumptions. 

In the opposite regime, R'g < Rg, the particles need 
to overcome a significant energetic barrier before they 
start feeling the attraction (see Figure The existence 
of this barrier is the major reason why we suggest to 
use linkers of at least several persistence lengths. In this 
case, the above Gaussian description may be applied. As 
we have shown, when the gyration radius of the linkers 
matches the scale of the repulsive potential, the particles 
may create a bound state without the need of overcoming 
any barrier. Thus, if the linkers are double-strand DNA 



molecules, the minimal value of 



lOOnm. 



The fundamental time scale of the problem is set by 
the lifetime of the AB contact, which can be estimated 



as 



vdR'' 



E 



(10) 



Here rj is the solvent viscosity, and we have assumed that 
the trial frequency is limited by the particle diffusion 
rather than by the departure rate of a linker-DNA. Since 
this desorption rate can be estimated as fcTc^^/ry, our 

regime requires > R'g~'^- This corresponds to a mod- 
erate marker-linker affinity, e/kT < 3\og{Rg/a) ~ 20 
(here a ~ 0.1 is the microscopic scale defined by 
the rigidity of hydrogen bonding between complementary 
DNA threads). 



Since the relaxation time (|l^) grows fast with the liner 
scales of the problem, the optimal regime corresponds 
to minimal value of Rg ~ Cmin — -ifJ-m, i.e d ~ Ifim. 
Within the assumption that the DNA-based key-lock 
complex would remain functional up to the coverage of 1 
marker molecule per ~ lOnm^, A^max ~ 10*. According 
to Eqs. (0-|), this yields ca/c^a = cb/c^b ~ 10~^ and 
cab/c^b ^ 10^^. The fundamental time scale ( p^ ) at 
the optimal regime is less than a minute. However, the 
true relaxation time of the system is determined by slow 
aggregation and growth processes. Hopefully, they both 
can be facilitated by using a commensurate substrate. 

Now that we have identified the optimal regime and in- 
terpreted the control parameters of the phase diagram in 
terms of the experimental ones, we should also determine 
how damaging is disorder for the expected behavior. The 
effect of the particle polydispcrsity is likely to be similar 
to the one in conventional colloidal self-assembly. The 
new effect is that the interactions are due to DNA chains, 
whose number per contact is of order of 1. Nevertheless, 
because of the reversibility of DNA-particle contact, the 
time-averaged interaction free energy varies only due to 
the discreetness of the markers. The corresponding dis- 
oerder can be estimated as SE/Uq ~ Rg^/Nmax/d < 0.1. 
As one can see from Figure |^, such a variation of E/Uq 
is tolerated in this problem, especially for the Gaussian 
potential. 
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